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The Comment by Markenscoff that criticizes a recent dynamic extension of the Peierls-Nabarro equation [Y- 
P. Pellegrini, Phys. Rev. B 81, 024101 (2010)] is refuted by means of simple examples that illustrate the interest 
of using an approach based on generalized functions to compute dynamic stress fields. 
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I. INTRODUCTION 

A Comment by Markenscoff 1 [hereafter referred to as (M)] 
criticizes several aspects of my paper 2 [referred to as (P)] 
on the derivation of dynamic extensions to the static Peierls- 
Nabarro equations. One remark 1 concerns my account of the 
author's work in the concluding Section of (P). Admittedly, 
the expression "extremely singular" I employ there is im- 
proper, for it might be interpreted as a statement that the prob- 
lem involves strong non-regularizable singularities, which it 
does not. To address the remaining points, I refer to Eshelby 's 
seminal paper 3 on dynamic dislocation motion by (E), and use 
equation numbers preceded by 'P' , 'M' of 'E' to refer to equa- 
tions in (P), (M), or (E). Unless otherwise stated, other equa- 
tion numbers relate to equations in the present Reply. 



n. DISCUSSION 
A. Eshelby's dynamic PN equation 

Ref. (P) focuses on obtaining dynamic equations of the 
Peierls-Nabarro (PN) type (to be solved numerically) to deter- 
mine the time-varying slip r)(x, t) on the slip plane y = of 
a planar dislocation with extended core, subjected to an arbi- 
trary time-dependent loading. The periodic pull-back force is 
arbitrary as well. Accordingly, no assumption on the slip func- 
tion must be made except smoothness. In this context I claim 
in (P) is that Eshelby's dynamic PN equation for screws (E21) 
misses the term (proportional to) —drj/dt featured by (P35). 
This is seen most directly from Appendix CI of (P), in which 
what I call (after Rosakis 4 ) "Weertman's equation", namely, 
the constant- velocity limit (P45) of the dynamic PN equation 
in the co-moving frame, is easily retrieved from (P35). This 
cannot be done with (E21), because the term —drj/dt proves 
essential in the process. 2 Whereas the kernel K(x,t) repre- 
sents the part of the cylindrical stress waves generated by the 
dislocation, that sweeps the slip plane, —drj/dt accounts for 
outward emission of updating waves collected at their incip- 
ience. On the basis of its negative sign 4 it can be interpreted 
(superficially, at least) as a radiative loss term. However, as 
shown below in Section HTT31 its operative role is to prevent a 
spurious stress term to show up at the dislocation core. 

Although claiming in (M) that (E21) is correct, Marken- 
scoff offers no convincing supporting argument. The self- 
force used by Al'shitz et al., namely (E28), is derived by Es- 



helby from (E26) as part of the equation of motion (EoM) that 
governs the time evolution of the dislocation position. This 
EoM assumes an arctangent slip function. Eshelby does not 
derive (E26) from (E21). Since (P) does not address EoM 
issues (dealt with elsewhere 5,6 ) but a more general question, 
the paper by Al'shitz et al. put forward in (M) is irrelevant 
here. It should be noted, however, that (E26) can be ob- 
tained from (P35). 6 Moreover in retrieving, 7 by assuming an 
arctangent slip function, the specific solution to Weertman's 
equation that corresponds to a sine pull-back force — rather 
than the latter general equation per se, contrary to what (M) 
states — Markenscoff appeals to non-zero y values. 7 This is 
not allowed in the context of (E21), so that this calculation is 
inconclusive with regard to the latter equation. 



B. On single-valuedness and distributions 

Prior to (P), DeWit 8 emphasized that the multivalued char- 
acter of static displacements associated to individual disloca- 
tions can be disambiguated by adding a distribution 9 (or "gen- 
eralized function" 10 ) part located on a Volterra cut. Whereas 
using a single- or a multi-valued displacement field can be 
seen in statics as a matter of choice, 11 made possible by an 
underlying gauge-invariant structure, 12 a single-valued treat- 
ment such as that used in (P) is perhaps more natural in dy- 
namics. Indeed, it is unclear how to extend to dynamics the 
physical arguments that have been put forward in statics to 
justify, in connection with mutivaluedness, the arbitrariness 
of this discontinuity surface. 12 Also, the surface spanned by a 
moving individual dislocation line, identified in (P) with the 
cut, becomes in dynamics a physical observable related to the 
knowledge of past trajectory. The latter is inherent to dynam- 
ics since past motion is used to compute retarded fields. 

Obviously, the necessary term —drj/dt is linked to this 
choice of jumping plane. The question of whether the cut 
might, after all, be arbitrary in a dynamic theory of disloca- 
tion lines should be examined carefully, perhaps in the light of 
gauge-related considerations^ 2 *^ Fig. 1 in (M) suggests such 
an arbitrariness, but it is noted that the result it relates to con- 
cerns the short-time response of a Volterra dislocation. At any 
rate, the choice made in (P) is appropriate 14 to extended planar 
(i.e., Somigliana) dislocations such as in the PN equation. 

Inasmuch as one adopts Mura's eigenstrains approach^ as 
in (P) the cut constitutes the support of the plastic eigenstrains, 
which must be properly accounted for when computing fields 
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on the slip plane. 16 Distribution theory can then be appealed 
to in order to make calculations without the need for a limiting 
process across the surface of discontinuity. 17 In (P) in partic- 
ular, Fourier transforms (FT) are always implicitly interpreted 
as generalized functions. 



C. Inclusion vs. dislocation formulations 

Two writings of the displacement u generated by an eigen- 
distortion (3*, representing an inclusion embedded in a host 
medium of infinite extent, are possible. To follow (M), we 
call them the "inclusion" and "dislocation" representations. In 
a plane problem these are, with Ax = x — x' and At = t — t', 



uf c (x,t) = -Idt'dyGij (Ax, At)C jklm pt mtk (x', 

uf s (x,t) = -[dt'd 2 x'G ij!k (AyL,At)C jklm f3f m (x',t')l2) 

in which Gij is the Green's function, Cy/./ is the tensor of 
elastic moduli, and where integration over x' is on infinite 
space. Their eventual reduction to "surface" integrals depends 
on the precise form of /3*. As far as x' is concerned and if f3* 
is bounded by some constant, t and t' stand as mere param- 
eters, and boundary contributions at infinity can be ignored. 
Equation (M8) is obtained from (O by Mura 15 (see p. 351). 
Expressions (fTJ and (fJJ differ only by a partial integration 
with respect to x', which makes them formally equivalent. A 
FT makes this obvious^ 

/f d 2 k 
dt'J ^-^fcfcGy (k, t~t')C jk i m /3i m (k, t'), 

(3) 

and the factor \k k can be interpreted as a derivative pertaining 
either to Gy or to /3f m . Calculations have been done in (P) us- 
ing ([3), but it is instructive to compute itj in direct space with 
(Q]l and (0. All formulations are found equivalent in the tran- 
sient regime contrary to what the first paragraph of (M) sug- 
gests. To show this, consider a time-varying eigendistortion 
localized on a surface S of normal n(x) by means of a sur- 
face distribution 5s(x), suitable to represent dislocations. 15 
Specifically, /?*j(x, t) = bi(x.,t)n,j(K)5s(x), in which &j is 
the local Burgers vector. Then, in the "inclusion" formulation 
the source term is the sum of one single- and one double-layer: 



d k p*j = d k {binj)8 s {-x) + hrij d k 8 s {x.). 



(4) 



For a Volterra screw of slip plane y = one has f3*j (x) — 
b(t)5a8 j28(— x)8(y) in which 9 is Heaviside's function. The 
part of <5g(x) is played by 8(y), and that of b by fr;(x, t) = 
b{t)9{-x)5 l3 . Thus 

S fc /3*.(x,t) = b(t)S i3 8 j2 [-5(x)5(y)5 kl + 9(-x)S'(y)S k2 ]. 

(5) 

For a screw in isotropic elasticity the single-layer contribu- 
tion vanishes (p, is the shear modulus): Cj k i m d k (3i m ('x, t) — 
/ib(t)9(—x)S' (y)Sj 3 . Focusing on the elementary time-depen- 
dent solution (P28), which lies at the root of the derivation of 



the term —dr\jdt in (P35), requires taking b(t) = bS(t). The 
distributional derivative in 5' (y) then entails 

u l 3 nc (x,t) =-nb [ dt'8{t') [ dx' [ dy' 



G 33 (x -x',y- y', t - t')5'(y') 



libdy f dx' G 33 (x - x ,y,t). 



(6) 



On the other hand, the "dislocation" formalism directly yields 



uf s (x, t) = -[ib I dx' d y G 33 (x - x', y, t). 

J — OO 



(7) 



Whereas in the "dislocation" formulation, by definition, the 
derivative d y takes place inside the integral over x', an am- 
biguity might subsist in the "inclusion" formulation, because 
whether the result should involve d y J dx' . . ., or J dx'd y . . ., 
or either indifferently, depends on whether the integrals over 
x' and y' in (|6]l can be harmlessly interchanged or not. In writ- 
ing the final expression of (|6]l it has been arbitrarily assumed 
that the integral over x' was done first. Simplifying the writ- 
ing of the integrals in the right-hand sides (rhs) of © and © 
by an obvious change of variables in x', and introducing 



G 33 (x,y,t) 



+ 00 



dx' G 33 (x',y,t), 



the proof of the equivalence consists in showing that 

f + OO 



d y G 33 (x,y,t) 



dx' d y G 33 (x',y,t), 



(8) 



(9) 



where, with r = \J x 2 + y 2 , the Green's function reads^ (cs 
is the shear wave velocity) 



G 33 (x,y,t) 



2^ (Cst ~ r )+ ■ 



(10) 



In this writing, the generalized function a;" = x a if x > 0, 
and otherwise, has been introduced to denote what can be 
written 9(4t 2 - r 2 )/^/c|i 2 - r 2 . Then, for c s t ^ \y\, 



G 33 {x,y,t) 
6(t)c s 



0(t)c s 

277/1 



27T/1 

0(t)c s 



0(c s t-\y\) 



dx (c s t — x — y 

a-\u\) 



2\-l/2 



du- 



(ii) 



tt9d{x, y, t) + 9(cst — r) arccos 



V4t 2 - v 2 



where 9d(x, y, t) stands for the characteristic function of the 
spatio-temporal domain D delimited by the simultaneous con- 
straints {x < 0, r > est, \y\ < est}, equal to 1 in D and to 
elsewhere. The function G 33 is represented in Fig. [T]at t = 1 
(with fj, = 1, cs = 1). The plateau is supported by domain 
D. On the circle x 2 + y 2 = c|t 2 , G 33 (x,y,t) is equal to 
9(— x)/A and connects continuously to the plateau on the left 
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mutually in the rhs. Owing to the delta-sequence 

-a - e 2 ) 



lim 



-ttS(o), 



(15) 



one finds lim f Re / 



-^(a 

a v 



-1/2 



Using this result in (fT2l with a — c\t 2 



2n9(-x)5(a). 
y 2 , shows that Ui 

in the dislocation formulation (0 again yields (P28). Thus, 
this alternative derivation produces the same result, which il- 
lustrates the consistency of the approach. 



FIG. 1. Function Ga(x,y,t) &tt — 1. 



D. The Dirac term in Equ. (P49) 



semi-circle (x < 0). Discontinuities are present along the 
Ox axis at x = for y = ztc^t, and along the Oy axis at 
y = ±cst for x < 0. The latter are responsible for the Dirac 
term S(cst — \y\) in (P28), this equation being retrieved from 
taking the y-derivative of the above expression of G33. The 
above shows that the Dirac sheets at y = ±c$t arise for x < 
as the envelope of circular waves emitted at each point of the 
perturbed region {x < 0,y = 0}. 

Turning now to the "dislocation" formulation, one has^ 



o+y 



ax. 



(—1 < a < 0) in the sense of general- 



ized functions, or more formally a Pf x°^ with Pf denoting 
Hadamard's finite part, in the sense of distributions (which re- 
quires using test functions). Then, for est ^ \y\, 



-I- 00 



dx' d y G 33 (x' \y,t) 



+00 



d ( f ) C S„. I j J , 2,2 /2 ,,2^-3/2 



V / dx (c s t - x -y 



(12) 



The integral in (fT2l can be computed conveniently without 
appealing to test functions, by using a representation of the 
generalized function similar to those found in Ref. [To| 

(e 2 t 2 -x 2 - y 2 y + 3/2 = lim Re[c s 2 t 2 - (x + ie) 2 - y 2 }^ 2 . 

(13) 

In the limit, the real part operator Re sets the rhs. to zero for 
Cgi 2 — x 2 — y 2 < 0, as must be. Carrying out the integral 
yields, for any nonzero real constant a, 



+00 



dx' 



1 



[a - {x> + ie) 2 ] 3 / 2 a 



x + le 



+ 0(-x)0(-a- e 2 )- 



2c 



y/a — (x + ie) 2 
(e > 0). (14) 



We remark that integral / is not defined for x < and 
a = —e 2 since in this case the integrand blows up as |x'|~ 3 / 2 
near the origin. Moreover, Re (x + ie) / Jo - (x + ie) 2 is, for 
a < — e 2 , a discontinuous function of x at x = because of 
the cut in the principal determination of the square root. Equa- 
tion ( TBI can be verified by taking derivatives with respect to 
x of both sides: then, the Dirac term created by the latter dis- 
continuity and that coming from the rightmost term cancel out 



Turning now to (P49), 19 Markenscoff's objections against 
the Dirac term <5(£ — x)£, where £(t) is the dislocation posi- 
tion, misunderstand its nature and working. In the formulation 
of (P49), mathematically correct in the sense of generalized 
functions although perhaps inconvenient for numerical calcu- 
lations, the Dirac term (which comes out of —drj/dt) com- 
pensates for the singularity as r — » t and makes the result 
finite in the subsonic regime. It serves the same purpose as 
the subtraction and addition of a compensating term in (M9). 
A notable difference is that it arises here as a direct conse- 
quence of the formalism employed. 2 Regularizing devices are 
not unique, and (M9) is by no means the sole manner of writ- 
ing the elastic strain. In (P49), the current instant r = t should 
be approached in the integral by a limiting process, for in- 
stance by multiplying the integrand by e~ £ / ( -*~ r - ) and letting 
e — » 0. This is related to the equal-time value of the Green's 
function G 3 3(x — x',t — t') being defined only as a limit 18 
t t' + (instead, considering the out-of-plane stress at small 
y 7^ would induce a natural time cut-off of the correspond- 
ing integral of order e ~ \y\/c$). With this precaution one can, 
e.g., deduce from (P49) the stress on the slip plane generated 
by a Volterra screw that at t = jumps from rest to a constant 
velocity v, a prototypical instationary case. The result reads 6 

(cli 2 -z 2 )f " 



2tt 



.Volt 



fib 



(x,t) = 



1 



1 



cs 



p.v. 



vt 



- 7rsign(w)[(t>/c s ) 2 - l]]/ 2 6(x - «t).(16) 

There is no Dirac term at the dislocation position in the sub- 
sonic regime \v\ < cs, as must be. The above-discussed com- 
pensation has occurred in the course of the derivation. 6 Again, 
the criticism in (M) is unjustified. Had the term —drj/dt not 
been present in (P35), and consequently no Dirac term fea- 
tured by (P49), a spurious stress contribution +tt(v/cs,)5(x — 
vt) would have remained in It should be noted that 

([T6l > embodies a former result for the subsonic regime by 
Markenscoff 20 [Equ. (17) in that reference] which it almost 
matches, and one for the supersonic regime by Callias and 
Markenscoff 21 [Equ. (2) at z = in that reference]. I write 
"almost" because the calculation 6 from (P49) to dT6l > automat- 
ically provides the principal value prescription p.v. in the first 
bracketed term of ([T6l >. which is absent from Markenscoff's 
subsonic expression. 20 This prescription regularizes the vicin- 
ity of the dislocation position x = vt. It makes the elemen- 
tary stress cr Volt a well-defined generalized function suitable to 
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convolution with a smooth core shape function (remark that f s (x) 
cr Volt is regular at x = for t > 0), as in statics (v — 0). 



Indeed, it is well-known that the Volterra stress kernel in the 
static PN equation is proportional to p. v. 1 /x, and not to 1/x 
(or 1/r) as stated in (M). 



E. The matter of Equ. (Mil) 



2tt, 
fib 



[a Volt */ e ]0M)=Re 



le 



v y 1 c\t 2 — (x + k) 5 
cs (x + is) — vt 



(17) 



By taking e a small number, dTvT > provides a numerical rep- 
resentation of (Tl~6b for any velocity. Moreover, the rhs. of 
(fTTI i. multiplied by — b/(2ir) (the minus sign being due to a 
different choice of dislocation sign), coincides with the sum 
of expression (74) in Ref. 7 and of — [b / (2ir)]x / (x 2 + e 2 ) 
in the subsonic regime for which this expression (74) holds. 
This sum of terms stands for lim^n [<9y M 3 0lt ('i Ui * fe] ( x )^ 
This illustrates the fact that in general, and contrary to what 
(Mil) states, the slip-plane limit can be interchanged with the 
convolution of a Volterra solution by a smooth core function, 
provided the Volterra solution is interpreted as a generalized 
function as in (P), in Ref. 6, and above. 



This leads us to equation (Mil). It originates from the 
strain of a Volterra screw moving at velocity v being pre- 
tended not integrable at x — vtl This cannot be, owing to the 
principal value prescription that should be there (see previous 
Section), and (Ml 1) should be an equality, rather than a non- 
equality, as physics commands (otherwise, even Weertman's 
stationary PN equation would be meaningless). To illustrate 
this point, a straightforward calculation shows that the con- 
volution of cr Volt as given by ( fT6b with the function (e > 0) 



III. CONCLUSION 

As far as explicit results are concerned, all of the par- 
ticular consequences worked out so far from the dynamic 
PN equations in (P) are in agreement with results obtained 
by other methods. Comment (M) ignores the generalized- 
function character of the Volterra solutions that legitimates the 
approach used. The criticisms against results in (P) have been 
convincingly refuted by means of simple explicit examples. 
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